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Quasi-Static Loading and Rupture 



Loading Phases: 

• 0) to A) - Quasi-static (QS) loading 

• A) to B) - Dynamic response 

Force 



Failure Criteria and Material Degradation 



Progressive Failure Analysis 


E 

Elastic 

property 



Benefits 

• Simplicity (no programming needed) 

• Convergence of equilibrium iterations 

Drawbacks 

• Mesh dependence 

• Dependence on load increment 

• Ad-hoc property degradation 

• Large strains can cause reloading 

• Errors due to improper load redistributions 
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Failure Criteria and Material Degradation 



Progressive Failure Analysis 
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Progressive Damage Analysis - Regularized Softening Laws 
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Strength-Dominated Failure 
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Before damage 


F = Act = EA 
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After damage 

F = Act = EA 


For stable fracture under A control: — 
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For “long” beams, the response is unstable , dynamic, and independent of Gc 
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Fracture-Dominated Failure 
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Crack propagates unstably once driving force G(a, a 0 ) reaches G lc 


6 


Fracture' Dominated Failure 
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Crack propagates stablv when driving force G(a, a 0 ) > G Init 
Unstable propagation initiates at G Init <G<G C 


Mechanics of Crack Arrest 





Crack arrest due to decreasing G 
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Mechanics of Crack Arrest 




Large strain rates often result in lower fracture toughness 
and delayed arrest 



Griffith Criterion and Stability 


Griffith growth criterion 





Stability of equilibrium propagation 
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< 0 unstable 


no growth 
equilibrium growth 
dynamic growth 


Wimmer & Pettermann 
J of Comp. Mater, 2009 
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Stability of Propagation with Multiple Crack Tips 




Scaling: The Effect of Structure Size on Strength 




Scaling from test coupon to structure 



Range of main practical 
interest for size effect 


Extrapolation that must he 
anchored in a theory 
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Cohesive Laws 



Bilinear Traction-Displacement Law 



Two material properties: 

• G c Fracture toughness 

• <J C Strength 


Yield or Strength Criteria 



log D 

Characteristic Length: 
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Crack Length and Process lone 


Force, F 



Brittle: 

a 0 > 100 l p 

Quasi-brittle: 

100 l p >a 0 > 5 l p 

Ductile: 



Applied displacement, A 


5 l p >a 0 
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Crack Length and Process lone 


Force, F 


LEFM error 




Applied displacement, A 


Brittle: 

a 0 > 100 l p 

Quasi-brittle: 

100 l p >a 0 > 51 p 

Ductile: 

5l p >a o 
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Strength and Process lone 



Force, F 


G. = constant 


Decreasing 


LEFM error 




As the strength a c decreases, 

1 . the length l of the process 
zone increases 

2. the error of the Linear 
LEFM solution increases 


Applied displacement, A 
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Size Effect and Material Softening Laws 



Damage Evolution Laws: 


Two material properties: 


Each damage mode has its 
own softening response 


• CJ c Strength 

• G c Fracture toughness 




log D 


Material length scale 
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Progressive Damage Analysis (Maimi/Camaniro 1007) 



Damage Modes: 
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LaRC04 Criteria 

• In-situ matrix strength prediction 

• Advanced fiber kinking criterion 
Prediction of angle of fracture (compression) 
Criteria used as activation functions within 
framework of continuum damage mechanics 
(CDM) 


Damage Evolution: 

Thermodynamically-consistent material 
degradation takes into account energy 
release rate and element size for each mode 


d, = 1 - j exp(/t,. (l - /. )) 

J i 

f{. LaRC04 failure criteria as activation functions 
i = F + \F~\M y+ \M y ~\M s 



Bazant Crack Band Theory: 
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Critical (maximum) finite 
element size: 
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log (Strength, MPa) 


Predicting Scale Effects with Continuum Damage Models 

<1 2mm tl 1mm il iUmu «l *mii» <1 lumm 

Prediction of size effects in notched composites 

• Stress-based criteria predict no size effect 

• CDM damage model predicts scale effects w/out calibration 

(P. Camanho, 2007) 

1 — 1 Hexcel IM7/8552 [90/0/45/-45] 3s CFRP laminate 





Deformed 
scale: 2x 




log (Strength, MPa) 


Process lone and Scale Effect in Open Hole Tension ^ 
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Stress distribution 

(P. Camanho, 2007) 
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Scale effect is due to 
relative size of process zone 
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Length of the Process Zone (Elastic Bulk Material) 



Short Tensile Test 

Lexan Polycarbonate 

Cohesive 
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Cohesive Laws - Prediction of Scale Effects 



• The use of cohesive laws to predict the 
fracture in complex stress fields is explored 

• The bulk material is modeled as either 
elastic or elastic-plastic 


Lexan Plexiglass tensile specimens (CT Sun) 



h/a=l (short process zone) 


l cz =4.7 mm. 

◄ ► 


Observations: 

• LEFM overpredicts tests for h/a<1 



h/a = 0.25 (long process zone) 


l cz = Width 
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Study of size effect: measuring the Recurve 


Double-notched compression specimens 
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Characterization of Through-Crack Cohesive Law 


Compact Tension (CT) Specimen 



Thin 

multidirectional 

laminate 






a 0 ' 


A a 


h w 




Experimental setup 



Characterization Procedure: 

1 . Measure R-curve from CT 
test 

2. Assuming a trilinear 
cohesive law, fit analytical 
R-curve to the measured V 
R-curve 

3. Obtain the cohesive law 
by differentiating the 
analytical R-curve 
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Size- Dependence of R-Curve 
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(a) Small CT 


(b) Large CT 
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Displacements 
measured through 
digital image 
correlation (DIC) 
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Plotting the R-curve as a function of 
the notch displacement removes the 
size-dependency 
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R-Curve Effect in Fiber Fracture 



Mode ll'Dominated Adhesive Fracture 
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Adhesive thickness: 0.13 mm 



Tip of adhesive 
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ENF J-Integral from DIG 





VIC3D output 
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Mode II cohesive law 




MMB Test - Analysis Results 
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Mixed mode bending (MMB) test fixture 


Nominally identical bonded MMB specimens sometimes fail in 
quasi-static mode and others dynamically. Why? 
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Double Delamination in MMB Tests 







Unexpected failure mechanism 

Two delamination fronts run in 
parallel: one in the adhesive, 
the other in the composite 


When the fiber bridge breaks, the crack grows unstably in the 
composite causing the drop in the load-displacement curve 
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Modeling the Double Delamination 



• A model was developed to evaluate the observed double 
delamination phenomenon 


• The model contains two additional cohesive layers within the 
composite arms 


MMB test specimen 
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Model of MMB specimen with double 
delamination 
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• This failure mechanism is often observed in bonded joints 
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Representative Volume Element and Micromechanics 


Why Micromechanics? 

Assumption: 

“Micromechanics has more built-in physics because it is closer to 
the scale at which fracture occurs ” 

Why NOT Micromechanics? (Representative Volume Element [RVE]) 

• Problem of localization 

• Randomness of unit cell configurations 

• Lengthscales missing 

• Characterization of material properties, especially the interface 

• Computational expense 
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RVE: 1) Problem of Localization 



Hardening 





Softening 







Localized 



Scale of RVE 
cannot be 
eliminated 
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RYE: 2) Randomness of Unit Cell Configurations 




Melro et al. IJSS, 2013. 


Fracture is a combination of interacting discrete and diffuse damage mechanisms 



Polymer Crack Bridging Plastically Deformed/ 
Matrix of Fibers Microcracked Matrix 


Bloodworth, V., PhD Dissertation, 
Imperial College, UK, 2008. 
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RVE: 3) Issue of Lengthi Scales 



RVE may not account for: 

• Ply thickness 

• Longitudinal crack length 

• Crack spacing 
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Transverse Strength of 90° Ply, GPa 


Matrix Cracking — In Situ Effect 




Potential crack plane. 



Thickness propagation 
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Transverse Matrix Cracks w/ One Element Per Ply 



Multi-element model: Conventional single-element: 

correct crack evolution no opening w/out delam. 



Modified single-element: 
correct Energy Release Rate 





Crack Initiation, Densification, and Saturation 



Van der Meer, F.P. & Davila, C., JCM, 2013 



Cohesive zone 
Traction-free cohesive zone 
Delamination 
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Material Inhomogeneity 



Initial crack density in a uniformly stressed laminate is 
strictly a function of material inhomogeneity 
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F Leone, 2015 
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Effect of Transverse Mesh Density on Crack Spacing 


F Leone, 2015 
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What Happened to Quadratic Convergence!!?? 


Commercial finite element vendors and 
developers are providing more and more 
tools for progressive damage analysis. 


But, if the load incrementation 
procedures do not converge... 


... more analysis tools 
more rope! 
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Techniques for Achieving Solution Convergence 


• Viscoelastic Stabilization 

• Delayed damage evolution 

• Implicit dynamics or Explicit solution 

• Arc-length techniques 

• Dissipation-based arc-length 
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Gutierrez, Comm Numer Meth Eng (2004) 
Verhoosel et al. Int J Numer Meth Eng (2009) 


Constant energy 
dissipation in each 
load increment 
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QS Solution of Unstable OHT Fracture 




Van der Meer, Eng Fract Mech, 2010 
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Open Questions 


* Is the QS solution physical? 

* Are the dynamic effects necessary? 

* Which solution provides more 
insight into failure modes? 




Implicit 


Explicit 
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Concluding Remarks 



• A typical structural tests usually consist of three stages: 

1 . QS elastic response without damage 

2. QS response with damage accumulation 

3. Dynamic collapse/rupture 

• Most structural failures exhibit size effects that depend on load 
redistribution that occurs during the QS phases 

• Correct softening laws based on strength and toughness considerations 
are required 

• Dynamic collapse/rupture is a result of the interaction between 
damage propagation and structural response 

• A stable equilibrium state often does not exist after failure under either 
load or displacement control 

• Onset of instability (failure) occurs when more elastic strain energy can 
be released by the structure than is necessary for damage propagation 

• Simulation of unstable rupture is often needed to ascertain mode of 
failure and to compare to test results 
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